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1 Introduction and main results 

Let us consider a standard 2n x 2n Gaussian Unitary Ensemble (GUE) of Hermitian random 
matrices W n : 

w n = wl (W n ) xy -L=(&* • /,/,,,). (1) 



2n 

where £ xy , r) xy , x,y = — n, ... ,—1,1, ... ,n are i.i.d. Gaussian random variables with zero mean 
and variance 1/2. Consider also the normalized eigenvalue counting measure (NCM) N n of the 
ensemble (jT)). defined for any Borel set A C R by the formula 

„„ (A ) = #%i^i, (2 ) 

In 

where Aj, i = 1, . . . , 2n are the eigenvalues of W n . 

Suppose now that ensemble (^Q) has also an additional symmetry of negative (positive) indices 
x and y. We consider four different cases of symmetry: 

1- (W n ) xy = (W n )_ y _ x , (3) 

2- (W n ) xy = (W n )_ x _ y , (4) 
3. (W n ) xy = (W n )_ xy , (5) 
4- (W n ) xy = {W n ) y _ x . (6) 

The Gaussian unitary ensemble and Gaussian orthogonal ensemble (GOE) was considered 
in numerous papers (see e.g. [3]). The Gaussian unitary ensemble with additional symmetry of 
type © was proposed in the papers 00] as an approach to the weak disorder regime in the 
Anderson model. This ensemble was also considered in the papers [2115]. In all these papers 
ensemble © was called as flip matrix model and studied by some supersymmetry approach and 
moments method. In this paper an approach is proposed that is simpler and the same for all 
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four cases This approach is a version of technique initially proposed in [7j and developed 

in the papers H EJ IE] ■ 

Using this technique we obtain the following results. 

First two ensembles (j3J) and @ are GOE-like. 

Proposition 1. The NCMs Nn^ and Nn^ of the ensembles © and (JI} converge weakly with 
probability 1 to the semi-circle law N sc 

N sc (d\) = (2^rV4-A 2 x h2 , 2 ](A)dA 
and the n~ l -asymptotics of the correlation functions 

FU( Zl ,z 2 ) = E { (gM( Zl ) - EgM(zij) (g%\z 2 ) - Eg${z 2 j) } , i = 1, 2 
of their Stieltjes transforms 

J-oo * — Z 

coincide with corresponding 2n x 2n-GOE asymptotic (2n)~ 2 S(zi, z 2 ) @].' 
F«(*i,z 2 ) = (2ny 2 S(z 1 ,z 2 ) + o(n~ 2 ) , 

q( 9 = ? / / ac (zi) -/ ac (z 2 A 2 / 7 v 

l l, 2j ; ' { ) 

W>0 



where 



is the Stieltjes transform of the semi-circle law N sc . 
The fourth ensemble © is GUE-like: 

Proposition 2. The NCM of the ensemble (JHJ) converges weakly with probability 1 to the 
semi-circle law N sc and the n _1 -asymptotic of the correlation function of its Stieltjes transform 
coincides with divided by 2 (i.e. GUE asymptotic). 

As for the third ensemble, the additional symmetry produces new limiting NCM and corre- 
lation function: 

Theorem 1. The NCM Nn of the ensemble (JSJ) converges weakly with probability 1 to the 
limiting non-random measure N 

N(dX) = \s(X)dX + i- V6 - (A 2 - A- 2 )x [ -A +! -A_]u[A_,A + ](A)dA ! (8) 



where \± = \/3± 2y/2 and the n 1 -asymptotic of the correlation function of its Stieltjes trans- 
form is given by the formula 

4" (a. a) = ( 2 «r 2 22) + o (>r 2 ) , 

oh, - (2 /<;•> ^ 2 + f i?' + ■> ) n +</<*»-'. w 

V - Z2) 2 2z(zif{z 1 )f(z 2 ) J k t^ 2 

where f(z) is Stieltjes transform of the limiting measure N . 

This result is somewhat unexpected for the Hermitian Gaussian random matrix ensemble 
with the rather large number (of the order re 2 ) of independent random parameters. But it shows 
how much the additional symmetry may affect the asymptotic behavior of the eigenvalues. 
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2 The limiting NCMs 

In this section we consider the limiting normalized countable measures of the ensembles (J3J)— (JBJ). 
In that follows we use the notations 



1 - 

G(z) = {W n - z)- 1 , g(z) = -Y, G J-^ 
1 n 



2n 

'j=-n 



9k/ 2n 

j=-n 

and (•) to denote the average over GUE. We also use the resolvent identity 

G(z) = -z- x I + z- x W n G{z) 

and the Novikov-Furutsu formula for the complex Gaussian random variable ( = £ + ir\ with 
zero mean and variance 1, and for the continuously differentiable function q(x,x) 

ECg(C,C) = Ej=g(C,C), (io) 

We will perform our calculations in parallel for all four ensembles. First, let us observe that 
properties ((5J) (JSJ) are valid not only for the matrices of ensembles © © but for their powers 
and hence also for their resolvents. Indeed, using induction by m and the symmetry of summing 
index we obtain: 

n n 

l=—n l=—n 
n 

= E ( w nh-j (wr 1 )^ = (wr 1 )^ . 

l=—n 

Thus, Gjk(z) = G^ k -j(z). 

n n 

2 - = E ( w nh (wr 1 )^ = E W;-* (^ m+1 u 

l=—n l=—n 
n 

= E (*c + v* = (wr 1 ).^, 

l=—n 

Thus, G jk (z) = G- j - k (z). 

n 

i=— n 
n 

Thus, 



i=— n 



G jk {z) = G_ jk {z) - z 1 (S jk - 5_ jk ) (11) 
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and, hence, 

n 

g(z) = g{z) - z" 1 , where g(z) = ^ G r _ 



\z . 



Unfortunately, there is no any such property for the fourth ensemble. 

Now using the resolvent identity for the average {G pq (z)), relation ()1U|) and formula for the 
derivative of the resolvent 

G'{z) ■ X = -G(z)XG(z), (12) 

we obtain 

(G pq (z)) = -z-% q + z- 1 {(W n G{z)) pq ) 

ft ft \ 

G r q{z f 



_! _! i ^ /i ( a , d \ 

~ z pq + z V^ r e r M\9C pr +z d Vpr ) 

= - z -H pq + z ^ (G rj {z)(W' n ) jk G kq {z)), 

r,j,k=—n 

where W' n = | — h ig^—^j W n . Now we calculate W' n for all four ensembles: 

j yyf _ 1 / 5jp5kr + Sjr^kp + Sj-r&k-p + $j-p5k-'. 

n 2\< r 2n V —djpdkr + Sj r 5kp — 5j- r 5k- p + 5j~ p 5k- 

{5jr$kp 4~ ^j—p^k—r) i 

V <zn 



\/2n 

( ^jp^kr 4~ fijrOkp 4~ fij—p^k—r 4~ ^j~r^k—p 

2\/2n V —bjpbkr + 5j r 5k p — 5j- p 5k- r + ^j-r^k-p 

— (5j r 5k p 4~ <5 j—r^k—p) j 



2 m/ _ 1 | fijpfikr 4" &jr&kp 4" 3j—p3kr 4" 5j r 5k—p \ 

n 2\/2n V —djpfikr 4- 5j>5fc P - 5j- p 5kr 4- 5j r 5k~ p J 

— /— — (.3jr$kp 4" Sjr^k—p) > 



^ ^ | $jp$kr 4~ &jr&kp 4~ $jr$k—p 4~ $j— p &kr 

n 2y/2n V —djpdkr 4- <Sj>5fcp - 5j r 5k- p + 5j- p 5k r 

(5j r Skp + Sj-pdkr) ■ 



>2n 

Using these formulas, we obtain the following relations: 



/ 1 n 

1. (Gpq(z)) = - z" 1 (g(z)G pq (z)) - z~ l { — ^ G r _ p (z)G_ r(? (z) 

r=— n 



2. (G M (*)) = -z" 1 ^ - z" 1 <^)G M (z)) - z" 1 (g(z)G. pq (z)) ; 

3. <G M (z)> = -z~H pq - z~ l (g(z)G pq (z)) - z- 1 (g(z)G- pq (z)} ; 

4. <G w (z)> = -z-H pq - z~ l (g(z)G pq (z)) - z~ l /i- £ G r _ p (z)G r(? (z)\ . (13) 

\ r=— n / 
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Now we put p = q in all four cases and p = —q another time in the second case, and apply 

n 

7^- ^2 . Thus, using also the additional symmetries of the resolvents of ensembles (JSJ)-©, we 
obtain: 



TrG 2 ( Z ) + <s°(z)sM) 



1. toM^-.-'ti + faM)')-.- 1 ^ 

where g°(z) = g{z) - (g{z)}; 

2. (g(z)) = -z~ l (l + (g(z)) 2 + (g 2 (z))) - z~ l [(g°(z)g(z)) + (g°(z)g(z))] , 
(g(z)) = -2Z- 1 (g(z)) (g(z)) - 2z^ (g°(z)g(z)) ; 

3. (g(z)) = -z- 1 (l + (g(z)) 2 + (g(z)) (g(z))) - z^ [(g°(z)g(z)) + (g°(z)g(z))] , 
g(z) = g(z) + z" 1 ; 

4. (g(z)) = -z' 1 (l + (g(z)) 2 ) - z~ l [i- (±. Tr P(z)G(z)^ + (g°(z)g(z)) 



(14) 



where matrix P(z) is defined by P xy (z) = G y - X (z). 

In the appendix we prove that the variances of random variables g{z) in all cases above are of 
the order 0(n~ 2 ) uniformly in z for some compact in C± (as well as the variance of g(z) in the 
second case). Besides, using Schwartz inequality for the matrix scalar product (A, B) = Tr AB, 
we obtain 



^Tr P(z)G(z) 
In 



< (i- Tr P(,)pt W ) V2 (_L MW Ot W ) < * 



2n 



Tr G 2 {z) 



< 



llm z 



2 ' 



Thus, all terms in square brackets in all four cases are at least of the order 0(n 1 ). Hence, in 
the first and in the fourth cases we obtain the following limiting equation: 



/(*) = (i + A*)) , 

which is the equation for f S c(z) — the Stieltjes transform of the semi-circle Law. 
Besides, since 



(15) 



9{z) 



1 



pq ~ 2n 



— Tr (W n G(z)) 
In 



< 



1/1 \ 1/2 



l/2\ i 1 V 1/2 

<(— Tr W n Wl) <1, 



then for all z with e.g. |Im 2;| > 3 uniformly in n we have in all cases 

\l + 2z~ 1 {g{z))\> 1 -. 
Thus, in the second case (g(z)} of the order (9(n~ 2 ): 

(g(z)) = -2Z- 1 (1 + 2z~ l (g(z)))- 1 (g°(z)g(z)) . 
Hence, the second case lead to the same limiting equation (fT5^l. 



(16) 
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As for the third case, it leads to the following equation 

f(z) = - z -i(l + 2f(z) + z- 1 f(z)). (17) 

Its solution in the class of Nevanlinna functions is the Stieltjes transform of the measure Q. 

The convergence with probability one in all four cases follows from the bounds for the vari- 
ances in the section bellow and the Borel-Cantelli lemma. 



3 The correlation functions 

As in the previous section, we perform our calculations in parallel for all four ensembles. 

Using the resolvent identity for the average (g (zi)G P q(z2)}, relations (fTU|) and 1[12|). we obtain 

1 n 

(g°(z 1 )G pq (z 2 ))=z 2 1 -= (9 (zi)Gr 3 (z 2 )(K) 3k G kq (z 2 : 

r,j,k=—n 
1 n 

+ z 2 3/2 (G lj (z l ){W^) 3k G kl {z 1 )G rq {z 2 )). 
y Zn ) l.r,j,k=-n 

Substituting in this relation the value of W n in all four cases and using the symmetries of the 
resolvents, we obtain 

1. (g°(z 1 )G pq (z 2 )} = -zz 1 {g { Zl )g{z 2 )G pq {z 2 )) - z^ 1 ( g° ( Zl ) ±-G 2 pq (z 2 ) 



i 1 

2 m \"2 



v 2nr V\ (G2(Zl)G(Z2) + (( G(Z2)G2(Zl) )^ P )) ; 

2. {g°( Zl )G pq (z 2 )) = -zz 1 (g°(z 1 )g(z 2 )G pq (z 2 )) - (g°(z 1 )g(z 2 )G- pq (z 2 )) 

~ H 1 -^ {((G\zi)G(z 2 )) pq ) + ((G(z 2 )G\ Zl ))_ p _ y )) ; 

3. (g°(z 1 )G pq (z 2 )) = -z? 1 (g°(zi)g(z 2 )G pq (z 2 )} - z? 1 (g° (zx)g(z 2 )G- pq (z 2 )) 

~ (((G 2 (^i)G(. 2 )) M ) + {{G\z x )G{z 2 ))_ pq )) ; 

/ 1 " 

4. (g ( Zl )G pq (z 2 )) = -z^ 1 (g ( Zl )g(z 2 )G pq (z 2 )) - z^ 1 ( 9°(zi) — £ G^ p (z 2 )G rq (z 2 

r=—n 



G\z 1 )G(z 2 )) p \ + (Y, {G\ Zl )) r _ p G rq {z 2 ) 



2 (2n) 2 

Then we put p = q in all four cases and p = —q another time in the second case, and apply 

n 

2^; X] an d obtain 

p=—n 

1. {g°{z l )g{z 2 )) = -2z^ (g(z 2 )) (g°(z 1 )g(z 2 )) 

~ *2 -1 7^2 ^TrG 2 ( Zl )G(z 2 )^ + n,„, (18) 

where 

n,n = ((g°(zi) { 9 °(z 2 )f) + i- ^ g °( Zl )±-TrG 2 (z 2 )^ ; (19) 
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2. {g { Zl )g{z 2 )) = -2z^ {{g{z 2 )) (g°(z 1 )g(z 2 )) + (g(z 2 )) (g°(z 1 )g(z 2 ))) 

i 2 / 1 



(g°(zi)g(z 2 )) = -2z 2 X {{g{z 2 )) (g°(z 1 )g(z 2 )) + (g(z 2 )) {g° ( Zl )g{z 2 ))) 



Z 



.! 2 / 1 



TVG 2 (z!)G(z 2 )) +r 3 , ; 



(2n) 2 \2n v y v 7 
where 

= -z 2 x ((g°(zi) (g°(z 2 )) 2 ) + {g°{z{) (g°(z 2 )) 2 )) , 

r 3>n = -2z^ 1 (g°(zi)g°(z 2 )g°(z 2 )) ; 

3. ^°(zi)^(z2)) = -4^ (5(^2)} <<7°(ziM* 2 )} - zf (g°(z 1 )g(z 2 )) 



-i 



1 / 1 



-2 



(2n) 2 \ 2n 



— Ti G 1 ( Z \)G(z 2 ) ) 



i ^(it(^)0(z 2 ))_J + 



r 4,' 

p=—n 



where 



r^ n = -2z 2 ~ 1 (g°{z 1 ){g {z 2 )) 2 )- 1 

4. (g°(z 1 )g(z 2 )) = -2z 2 x (g(z 2 )) (g° ( Zl )g(z 2 )) 



i 1 / 1 



{ 2nf\Yn^ G ^ G ^) + ^ 
where 

r 5 ,n = ~z 2 1 ^ (g°(zi)^ TrP(z 2 )G(z 2 )^ - z^ (g°( Zl ) {g°{z 2 )y 



Z7 1 



As we show in the appendix, all r Jjn , j = 1, . . . , 5 are of the order o(n~ 2 ). Thus, as one can 
easily show, all correlation functions F(z\,z 2 ) = (g° (zi)g(z 2 )) above are of the order 0{n~ 2 ). 
Moreover, since {g(z 2 )) is of the order 0(n~ 2 ) in the second case, its easy to see that cases one 
and two lead to the same relation for F{z\, z 2 ) 

F( Zl ,z 2 ) = -2Z" 1 (g(z 2 )) F( Zl ,z 2 ) - z 2 x -^ Q^Tr G 2 (z l )G{z 2 ) X j + o (n~ 2 ) . (21) 
As to the case four, it leads to 

F{z u z 2 ) = -2z 2 x (g(z 2 )) F(z u z 2 ) - ^"A; Tr G 2 ( Zl )G(z 2 )^ + o {n~ 2 ) . (22) 
Besides, due to the resolvent identity we have 

■^-T?G 2 (zi)G(z 2 ) = -J- (^TrG^) - ~ 9^) \ . (23) 
2n z\ — z 2 \2n z\ — z 2 ) 
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In addition, as we show in the appendix, in these cases 

i^'f (24) 

Thus, substituting in the relations (|21|). ((2~2*)) the expressions (|2*3)) . (|2*i)) and using the equa- 
tion l)15|) for the limit of (g(z)), we obtain in the cases one and two the GOE correlator asymp- 
totic and in the case four the twice less GUE asymptotic. 
To treat the third case we use and obtain that 

1 71 1 1 

- Y. (C 2 ( 2l )C fe ))_„ = -T rf - (zi )C( 22 ) + 

p=—n 1 

This gives the following relation for F(z±, z-i) 

- £L ( ±- + - J- (J- Tr G\ Z1 ) - ^-^ )) + („->) . ( 2 5) 
(2n) 2 zi-Z2\2n z x - z 2 // y 

We show also in the appendix that in this case 

1 rp ^2/ A (gO)) 1 ~ Z~ 2 , _ 2 , 

— Tr G (z) ) = , . .. 1 k + o in . 

2n w / z 1 + (g(z)) z~ l + z~ 2 v ; 

Substituting this relation in (|25j) we obtain 



F(z 1 ,z 2 ) = 



1 , 2 f(z!)-f{z 2 ) 
1 / (212 2 ) 2 *1 (zi-^2) 2 



?; 2 



1 + zf 2 + fcf'/fa) 

-4^7? n (i+^+4^))+°(«- 2 ). 

zi - z 2 ^1^2 k t^ 2 V ^fc / y 

Then, using the equation (|17|). we rewrite this relation in the form @. 



4 Conclusion 

The purpose of this paper was to answer the question: "Can the additional symmetry prop- 
erties influence on the asymptotic behavior of eigenvalue distribution of GUE?" The negative 
answer for the three cases of additional symmetry is not surprising, as these symmetries leave 
the number of independent random parameters of the order n 2 . The effect when in one case the 
additional symmetry essentially changes the limiting eigenvalue counting measure is very unex- 
pected, especially the appearance of the gap in the support of limiting NCM. Unfortunately, the 
physical application of this effect is unknown to the author, though one of the other considered 
ensembles (flip matrix model) was used as an approach to weak coupling regime of the Anderson 
model. 



A Appendix 

Proposition 3. The variance v = (^\g°(z)\' 2 ^ is of the order 0{n~ 2 ) in all four cases, and the 
terms rj >n , j = 1, . . . , 5 are of the order o(n~ 2 ). 
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Proof. First we proof that the variance is of the order 0(n 2 ) in all four cases. Indeed, in the 
first case, using with z 2 = z~[ = z, we obtain 

v(l + 2z' 1 (g(z))) = -z 2 1 -^ /-!- Tr G 2 (z l )G(z 2 )) + n,„. 



(2nY \2n 

Besides, using the Schwartz inequality we obtain from (fTT))) 
1 f _J_„_l - ^ 



ri,n < \Z\ tz -V H — 

VI Im ^l 2n\lmz\ 
Thus, due to the bounds lfT5)l and 



^TtG 2 ( Zi )G(z 2 ) 
2n 



< TT^' ( 26 ) 
1m z\ 



we have for |Imz| > 3 the inequality 

v < — + ^-v l/2 , 
i(2n) 2 2n 



which leads to v = 0(n 2 ). For the other cases the proofs are analogous. 
To prove r% n = o(n~ 2 ) for 



n, 



-Z2 1 (( 9 °(zi) (g°(z 2 )) 2 ) + i- (g°( Zl )±-TrG 2 {z 2 )^ 



we rewrite the second term in the parentheses as 

l / n , s l m ^, A l / „, , l „ 5 _ A 



Since the value (g°(zi)g°(z 2 )) is analytical and uniformly in n bounded for |Im 2:1,2 1 > 3, and 
since, due to the Schwartz inequality \(g°(zi)g°(z 2 ))\ < v = 0(n~ 2 ), its derivative on z 2 is also 
of the order 0(n~ 2 ) and hence the second term is of the order 0(n~ 3 ). 
To prove that the first term is o(n~ 2 ) let us consider 

(\9°(z)\ 4 ) = (i y g {z l )g°{z 2 )) 2 ) = (R°g(z 2 )} , R = (g°( Zl )) 2 g°(z 2 ), z x =~ 2 = z. 

Then, using the resolvent identity for the average (R°G pq (z 2 )}, relations (fTU|) and (|12|). we obtain 

1 n 

{R°G pq {z 2 )) = z 2 : 1 — (R°G rj {z 2 ){W' n ) ]k G kq {z 2 )) 

r,j,k=~n 

2 n 1 

+ *?TT\ W2 E (9 (z 1 )g (z 2 )G lj (z 1 ) {W^ jk G kl ( Zl )G rq (z 2 ; 

V 2n ) l.r,j,k=-n 

+ Z 2 1 -^sj 2 E ((^(^i)) 2 ^^) «) jfc c fei (^)G r(? (z 2 ; 

(2 n ) l.r,j,k=-n 

Substituting in this relation the value of WL and using the symmetry of the resolvent we obtain 



(R°G pq (z 2 )) = -z 2 x (R°g(z 2 )G pq (z 2 )) - z? ^R° ±-G 2 pq (z 2 )^ 
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-z 2 1 -^((g (z 1 )) 2 G s pq (z 2 )). 
Then we put p = q in all four cases and p = —q another time in the second case, and apply 

n 

2^; X] an d obtain 



p=—n 



f{z)\ 4 ) = (R°g(z 2 )) = -z^ 1 (R°g 2 (z 2 )) - z^i- (#1. Tr G 2 (z 2 ] 



^^n) 2 ( 9 ° {Zl)9 ° {Z2)l L TT ( G2 ^ G ^)) 



Using this relation, the bounds (|26|) and 

|4 



[R°g 2 (z 2 ))\ = \(R°g°(z 2 )g(z 2 )) + (R°g(z 2 )) (g(z 2 ))\ < , 

R°-L T rG 2 (z 2 )) <-%, 
2n I |Imz| 

we obtain that for \Imz\ > 3 (||g (z)| 4 ^ = 0(n~ 3 ). Thus, due to the Schwartz inequality the 

term (g°(zi) (g°(z 2 )) 2 ^ is of the order 0(n -5 / 2 ) and, hence, ri >n is of the same order. 

The cases of the terms Tj n , j — 2, . . . , 5 can be treated analogously, with exception for the 
last term of r§ n . The last term of (JUJ) 

can be treated as follows. 

First, observe that in the case four (g(z)) = o(n~ 1 ). Indeed, using (|TH|) with q = —p, we 
obtain 

(g(z)) = -z- 1 (g(z)) (g(z)) - z' 1 (g°(z) 9 (z)) - z~ l i- ^i-Tr {G(z)G T (z)) ^ , 

where G T is transpose of G. Due to the the Schwartz inequality for the trace, the last term in 
r.h.s. of this relation is of the order 0(n _1 ). Since the variance of g(z) is of the order 0(n~ 2 ), 
the second term is at least of the order 0(n~ 1 ) (in fact it is of the order 0(n~ 2 ), since, as one 
can show, the variance of g(z) is of the same order). Thus, {g(z)) is of the order 0(n~ l ). Its 
easy to show in the same way that 



3=-n 



is also of the order 0(n and its variance is of the order 0(n 2 ). 
Now, using the resolvent identity for the average of 



Y n E (G 2 { Zl )) p _ q G pq {z 2 ), 
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relations (fTU|) and 1)12(1 . we obtain 



<*> = "*2 1 (H*)) - ^7^2 ( E (G 2 (z 1 )) p ^ q G rj (z 2 ) (W' n ) jk G kq {z 2 ) 

^ / \p,q,r,j,k=— n I 



Z, 



P,q,r,3,K= 

1 1 / n 

q,r ,j,k,m=— n 



-^3/2\ E Gpj(zi){Wn) jk G(z 1 )kmG m - q (z 1 )G rq (z 2 ) 

1 



z, 

(2n) 



- — ^3/2 \ E C pm (zi)G(zi) mi (lO ifc G fc _q(zi)G r g(z 2 ) 

v / \p,q,r,j,k,m=—n 



Substituting in this relation the value of W^, we obtain 



(cD) = -z 2 l (h( Zl )) - z^ 1 (g(z 2 )<S>) - z -^-l— £ (G\ Zl )) p _ q G r ^ p (z 2 )G rq (z 2 ) 

\ p,q,r=—n 

2 1 i(i E (^i)g(4,gV,(^i) 

\ p,q=—n I 
\ p,q=—n I 



Z, 



z, 



The first term of the r.h.s. is of the order 0(n~ 1 ), the last five terms are of the same order, 
because of the Schwartz inequality and of the bounds for the variances of g(zi) and h(z\). The 
second term we rewrite as follows 

-z^ (g(z 2 )t>) = -z 2 l (g(z 2 )) <$) - z^ (g (z 2 )^) , 

where due to the Schwartz inequality the last term is also of the order 0(n~ 1 ). Thus, we conclude 
that (<5) is of the order 0(n _1 ) and, hence, the last term of r^ tTl is of the order 0(n -3 ). ■ 

Proposition 4. In the third case we have 

~ Tr G\z)\ = - {9{Z)) ? + o (n" 2 ) . 

2n V V z 1 + (g(z)) z- 1 + z- 2 V ; 

Proof. Indeed, we have 

Hence, we can just take the derivative of the identity l)14|) for {g(z)) . Since all terms of the order 
0(n~ 2 ) in square brackets in (|14[) are analytical and uniformly bounded for |Imzi2| > 3, they 
remain of the same order. Thus, we obtain the relation needed. ■ 
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